Breakdown of The Excess Entropy Scaling for the Systems with Thermodynamic 
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This articles presents a simulation study of the applicability of the Rosenfeld entropy scaling to 
the systems which can not be approximated by effective hard spheres. Three systems are studied: 
Herzian spheres, Gauss Core Model and soft repulsive shoulder potential. These systems demon- 
strate the diffusion anomalies at low temperatures: the diffusion increases with increasing density 
or pressure. It is shown that for the first two systems which belong to the class of bounded poten- 
tials the Rosenfeld scaling formula is valid only in the infinite temperature limit where there are 
no anomalies. For the soft repulsive shoulder the scaling formula is valid already at sufficiently low 
temperatures, however, out of the anomaly range. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 



I. INTRODUCTION 

It is well known that some liquids (for example, water, 
silica, silicon, carbon, and phosphorus) show anomalous 
behavior in the vicinity of their freezing lines jll-Eil. The 
water phase diagrams have regions where a thermal ex- 
pansion coefficient is negative (density anomaly), a self- 
diffusivity increases upon pressuring (diffusion anomaly) , 
and the structural order of the system decreases upon 
compression (structural anomaly) 6, 2]- The regions 
where these anomalies take place form nested domains 
in the density-temperature [S] (or pressure-temperature 
0]) planes: the density anomaly region is inside the dif- 
fusion anomaly domain, and both of these anomalous 
regions are inside the broader structurally anomalous re- 
gion. It is natural to relate this kind of behavior with the 
orientational anisotropy of the potentials, however, there 
are a number of studies which demonstrate the water- 
like anomalies in fluids that interact through spherically 
symmetric potentials 13" 31, 33-45|. 



It was shown IT], ll8[ that the thermodynamic and ki- 
netic anomalies may be linked through excess entropy. In 
particular, in Refs. 17, H the authors propose that en- 
tropy scaling relations developed by Rosenfeld 0, |43| 
can be used to describe the the regions of diffusivity 
anomaly. 

Rosenfeld based his arguments on the approximations 
of liquid by an effective hard spheres system. In this 
approach the kinetic coefhcients are expressed in re- 
duced units based on the mean length related to den- 
sity of the system d = p~^^'^ and thermal velocity Vth — 
{ksT /niY^'^ . The reduced diffusion coefhcient D* , vis- 
cosity Tj* and thermal conductivity k* are written in the 
form 



D* = D 



(1) 



,-2/3 



(2) 



(3) 



Rosenfeld suggested that the reduced transport coef- 
ficients can be connected to the excess entropy of the 
system Sex = {S — Sid) /{NkB) through the formula 



X = ax ■ e 



(4) 



where X is the transport coefficient, and ax and bx are 
constants which depend on the studying property (47| . 
Interestingly the coefficients a and b show extremely weak 
dependence on the material and can be considered as 
universal. 

Another expression for relating diffusion coefficient to 
the excess entropy was suggested by Dzugutov [4§]. In 
this approach the natural parameters of the system were 
chosen to be the particle diameter a and the Enskog col- 
lision frequency Fb = 4:a^g{a)p^/ ^^^'^ , where g{a) is the 
value of radial distribution function at contact. In case 
of continuous potentials the value of a corresponds to 
the distance of the first maximum of radial distribution 
function. Defining the reduced diffusion coefficient as 



n* — 
for it 



Dzugutov suggested the following formula 



n* 



0.049e^=, (5) 

where S2 is a pair contribution to the excess entropy 
(49I . It was shown that this relation holds for many sim- 
ple liquids. At the same time this equation is not strictly 
valid for liquid metals. In the work [50|] it was shown that 
in this case it is necessary to replace the pair entropy S2 
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by the full excess entropy Sex- It was also shown that 
Dzugutov formula does not work for silica modelled with 
an angular dependent potential [HO] ■ It allows to say that 
Rosenfeld relation is more general. 

Remember that original idea underlying the Rosen- 
feld relation is to refer the system under investigation 
to the hard spheres system. In this respect it is inter- 
esting whether the Rosenfeld scaling relation is valid for 
the systems essentially different from the hard spheres. 
One of the examples of such systems is the system with 
potentials with negative curvature 1^, 23] ■ It was shown 
in many publications that the behavior of such systems 
is very complex [4l|, 0, 0, [5ll - [57| . In particular such 
systems can form complicated structures, like cluster liq- 
uids or different crystal phases. They can demonstrate 
maximum on the melting line and reentrant melting and 
many other unusual properties. In particular systems 
with negative curvature can demonstrate anomalous be- 
havior 0, ii, [SSi, [H, m, [5^ . 

It was suggested that the Rosenfeld relations can hold 
even in the case of anomalous diffusion [l^, [ll]. For 
example, in the paper (l8| the dependence of both ex- 
cess entropy and diffusion coefficient on density are re- 
ported for the core-softening potential that consists of a 
combination of a Lennard- Jones potential plus a Gaus- 
sian well. This potential can represent a whole family 
of two length scales intermolecular interactions, from a 
deep double-well potential to a repulsive shoulder. Ac- 
cordingly to these dependencies both excess entropy and 
diffusion have non monotonic behavior which allows to 
preserve exponential dependence of the diffusion coef- 
ficient on the excess entropy. It means that the ther- 
modynamically anomalous regions are characterized by 
anomalous behavior of the excess entropy which induces 
anomalous diffusion as well. 

Another example of systems which can not be ap- 
proximated by a hard sphere model is the systems with 
bounded potentials j59H61l. I67i |. Since these potentials 
have no singularity in the origin the behavior of such 
system is strongly different from the behavior of hard 
spheres. 

One of the most common model with bounded poten- 
tial is the Gaussian Core Model (GCM). This system is 
defined by the potential 



Uair) 



ee 



(6) 



This potential was introduced by Stillinger [53] for sim- 
ulation of the plastic crystals system . The phase diagram 
of the GCM demonstrates two crystal phases - fee and 
bcc [i^]. Starting from the densities around pa^ « 0.25 
the melting curve has a negative slope. It was also shown 
that GCM demonstrates liquid state anomalies: density 
anomaly [H, [63[^_^ diffusion anomaly [gl, 65] and struc- 
tural anomaly '64]. Interestingly Stockes-Einstein rela- 
tion is also violated in the GCM system ^66i] . 



In the article [67[ an extensive study of another model 
with bounded potential was reported. This work is con- 
cerned to the Herzian spheres system which is defined by 
the interparticle potential of the form 



$(r) = 



0, 



r < a 
r > a 



(7) 



The phase diagram of the Herzian spheres system 
demonstrates very complex behavior, including many 
crystal phases and reentrant melting. Anomalous dif- 
fusion is also reported |67| . 

Taking into account that the behavior of the systems 
with negative curvature potentials and the systems with 
bounded potentials is rather different from the behavior 
of hard spheres a question arises if the Rosenfeld scaling 
relations are applicable for such systems. 

The purpose of this article is to analyze the validity 
of the entropy scaling for the systems with anomalous 
behavior. For our analysis we have chosen the diffusion 
coefficient since it is the simplest transport coefficient 
to calculate in simulation. The article is organized as 
follows. In the second section we describe the models in- 
vestigated in the present work and the simulation setup. 
Section III gives the results of the simulations and the 
discussion of these results. Finally the section IV repre- 
sents our conclusions. 



II. THE SYSTEMS AND METHODS 

Three systems were studied in the present work: 
Herzian spheres. Gauss core model and a soft repulsive 
shoulder system. 

For the simulation of the Herzian spheres we used 
a system of 1000 particles in a cubic simulation box. 
NVE MD simulation was carried. Equations of mo- 
tion were integrated by velocity Verlet algorithm. The 
time step was set to dt — 0.0005. The equilibration run 
was 5 • 10^ time steps and the production run 1.5 • 10^ 
time steps. During the equilibration the velocities were 
rescaled to keep the temperature constant. The diffu- 
sions were computed via the Einstein relation for the 
densities from p = 1.0 till p = 15.0 with the step 
Ap — 0.5. Additional simulations were done for com- 
puting the equation of state for the densities less then 
unity. Free energy of the liquid was calculated by inte- 
grating the pressure along an isotherm [s^] and the excess 
entropy was obtained from the relation S^^ — ^"^^^ , 
where U is the internal potential energy of the liquid. 
The simulations were done for the set of ten isotherms: 
T = 0.01; 0.02; 0.03; 0.05; 0.1; 0.15; 0.2; 0.25; 0.3; 0.5. 

In the case of GCM the system consisted of 2000 
particles. The time step was set to dt = 0.05. 
The equilibration and production runs were 4 • 10^ 
and 1 • 10® time steps respectively. The diffusion 
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was computed for the densities from p = 0.1 to 
p = 1.0 with the step of 0.1 for the isotherms 
T = 0.04; 0.05; 0.06; 0.07; 0.1; 0.2; 0.3; 0.4; 0.5; 1.0; 2.0. 
The calculation of the diffusion coefhcient and excess en- 
tropy was carried in the same way as for Herzian spheres. 

The last system considered in the present work is the 
continuous repulsive shoulder system introduced in the 
article 



41| |. The potential of this system has the form 



[1 - tanh{ko{r - (Ti))], (8) 



where /cq = 10.0. As it was reported in the paper [45| 
this system demonstrates anomalous behavior due to its 
quasibinary nature. Here we extend the investigation of 
the diffusion anomaly in the repulsive shoulder system 
with (Ti = 1.35 and check the Rosenfeld relations for this 
system in the anomalous region. 

For the simulation of the repulsive shoulder system we 
used parallel tempering technic [6^. The details of the 
simulation were described in [4^. We computed the dif- 
fusion coefficients along different isochors starting from 
the density p = 0.3 to p = 0.8 with the step 0.05. A 
set of 24 temperatures between T = 0.2 and T = 0.5 
was simulated. Taking into account the exchange of the 
temperatures at the same density more then a hundred 
runs at the same isochor was done. This allowed us to 
collect a good statistics on the temperature dependence 
of the diffusion coefficient. The diffusion coefficient along 
an isochor was approximated by a 9— th order polynome 
of the temperature. The excess entropy was calculated 
in the way described above. 

Usually excess entropy can be well approximated by 



the pair contribution only: Sex — Spair - 
where 



'S. 



s 



Spa^r = -\p J dr[g (v) Hg (v)) - (gir) - 1)], 



pair 1 



(9) 



where p is the density of the system and g(r) is the 
radial distribution function. We did not use the pair con- 
tribution to the excess entropy for the GCM and Herzian 
spheres because of the considerable overlap of the parti- 
cles for the bounded potentials. 

Since the potentials studied in the present work have 
negative curvature regions or are bounded they can not 
be approximated by an one component hard spheres sys- 
tem. It allows us to pose a question about the appli- 
cability of the entropy scaling to these systems both in 
Rosenfeld and Dzugutov forms. Note that Dzugutov re- 
lation (Eq. (5)) involves the size of the particles a which 
is ill defined for the negative curvature and bounded po- 
tentials systems. This makes problematic to apply the 
Dzugutov scaling rule to them. Because of this only 
Rosenfeld relations were used in this work. 



In this paper we use the dimensionless quantities: f = 
r/cr, P = Pa^/e, V = V/Na^ = 1/p, f = ksT/e. Since 
we use only these reduced units we omit the tilde marks. 



III. RESULTS AND DISCUSSION 

This section reports the simulation results for the dif- 
fusion coefficient and excess entropy of the three models 
described above and checks the validity of the Rosenfeld 
relation for these systems. 



Herzian spheres 

Low temperature behavior of the diffusion coefficient of 
Herzian spheres system was already reported in the work 
[67| . As it is seen from this publication the diffusivity 
shows even two anomalous regions at the temperature 
T = 0.01 where diffusion coefficient grows with growing 
density. In the present work the dependence of the dif- 
fusion coefficient on density along several isotherms was 
monitored. The simulation data are presented on the 
Fig. 1 (a) - (b). 

One can see from these figures that at low tempera- 
tures (Fig. 1(a)) the diffusion is non monotonic, while at 
high temperatures it monotonically decays with increas- 
ing density (Fig. 1(b)) and comes to a constant value 
(see inset of the Fig. 1(b)). 

It is worth to note that the meltin g te mperatures of 
Herzian spheres reported in the work 67| are of the or- 
der of 10~^, so the temperatures about 0.1 are extremely 
high for this model. This is easily seen from the Fig. 2 (a) 
- (b) where the radial distribution functions for the den- 
sity p = 6.0 are shown for the same set of temperatures. 
One can see that at T = 0.01 the liquid has short range 
structure which rapidly decays with increasing tempera- 
ture. At the temperature T = 0.1 the liquid looks almost 
like an ideal gas since g{r) comes to unity very quickly. 

Excess entropy also shows non monotonic dependence 
on density along an isotherm (Fig. 3(a) -(b)). One can 
see from these figures that at low temperatures excess 
entropy has two minima and a maxima in the investi- 
gated density range while at high temperature the first 
minima is depressed and the curves just change the slope 
smoothly. 

Now we turn to the Rosenfeld relation for the Herzian 
spheres system. The dependence of the reduced diffu- 
sion (see formula (1)) on the excess entropy along some 
isotherms is shown in the Fig. 4 (a) - (b). Looking at 
the curve for T — 0.01 (Fig. 4 (a)) one can divide it into 
three distinct regions with different slopes which we de- 
note as '1', '2' and '3'. The density increase corresponds 
to moving along the curves from right to left, i.e. region 2 
corresponds to the higher densities then 1, and 3 - higher 
densities then 2. As is seen from the plots the region 3 
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FIG. 1: DifTusion coefficient of Herzian splieres for a set of FIG. 2: Radial distribution functions of Herzian spheres at 

isotherms, (a) - T = 0.01; 0.02 and 0.03; (b) - T = 0.1; 0.2 p = 6.0 and a set of temperatures, (a) - T = 0.01; 0.02 and 

and 0.3. The inset of the figure (b) showes the high density 0.03; (b) - T — 0.1; 0.2 and 0.3. 
behavior of diffusion coefficient. 



rapidly disappears with increasing the temperature. Al- 
ready at T = 0.03 this region is negligibly small. Recall 
from the figures 1 and 3 that at low temperatures both 
diffusion and excess entropy behave non monotonically 
while at growing the temperature this effect disappears. 
This leads to the depression of the region 3 in the Fig. 4 
(a). 

The region 2 is rather stable. As one can see from the 
Fig. 4 (b) this region also becomes less developed with 
increasing the temperature, but it still preserves even for 
high temperatures. It makes the excess entropy scaling 



curve consisting from two parts of different slope and a 
cross region. 

Fig. 5 summarizes all the results obtained for Herzian 
spheres system. Ten different isotherms are shown there. 
As one can see only at the temperature as high as 0.5 the 
Rosenfeld linear relation between the logarithm of the 
reduced density and the excess entropy becomes valid. 
Remember that the melting temperature is of the order 
of 0.005, that is 100 times smaller. It allows to say that 
the Rosenfeld relation for diffusion coefficient of Herzian 
spheres is valid only in the infinitely high temperature 
limit. 
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FIG. 3: Excess entropy along and a set of isotherms, (a) 
T = 0.01; 0.02 and 0.03; (b) - T = 0.1; 0.2 and 0.3. 



FIG. 4: The dependence of reduced diffusion on the excess 
entropy (see formulas (1) and (4)) along some isotherms (a) 
- T = 0.01; 0.02 and 0.03; (b) - T = 0.1; 0.2 and 0.3. 



Gaussian Core Model 

The results for the GCM are qualitatively similar to 
the case of Herzian spheres. Because of this we do not 
explain them in detail. Fig. 6(a) and (b) shows the dif- 
fusion coefficient for the GCM system for a set of six 
isotherms. One can see that for low temperatures start- 
ing from the densities approximately 0.3 the diffusion 
coefScient demonstrate anomalous growth with increase 
of the density. We expect that at higher densities the 
curve bends downward but in the present study we have 
not measured so high densities for this model. At high 
temperatures the diffusion monotonically decreases with 



increasing the density. 

The structure of the liquid rapidly decays with the 
temperature increase. This is shown in the Fig. 7(a) - 
(b). One can see that at the temperature T — 1.0 g{r) is 
equal to unity almost in the whole range of the distances 
r. 

Like the diffusion coefficient the excess entropy has 
a minimum at the low temperatures, but with increas- 
ing the temperature it becomes monotonically decreasing 
function of the density (Fig. 8(a) - (b)). 

Fig. 9(a) - (b) demonstrates the Rosenfeld relation for 
the GCM for the same set of isotherms. Comparing to 
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FIG. 5: Rosenfeld excess entropy scaling of difFusion coeffi- 
cient for a set of ten isotherms for Herzian splieres system. 
The inset shows several high temperature isotherms in the 
enlarged scale. 



the case of Herzian spheres there is no region '3' in the 
low temperature curves of the GCM. We suppose that 
this region corresponds to the higher densities which we 
do not consider in the present study. At low tempera- 
tures most of the points belong to the region '2'. How- 
ever this region depresses with the temperature increase. 
Although even at the temperature as high as T = 2.0 
which is around 200 times higher then melting tempera- 
ture of GCM the curve still demonstrates a bend from a 
straight line at high densities. 

Finally Fig. 10 shows the whole set of isotherms in- 
vestigated in the present study. As one can see from this 
figure the system comes more close to the Rosenfeld re- 
lation with increasing the temperature. One can expect 
that at infinitely high temperatures the excess entropy 
relation is valid for the GCM, however in the tempera- 
tures range studied here, i.e. to approximately 200 -Tmeis 
a deviation from the linear behavior is still observed. 



Soft Repulsive Shoulder Model 

The last model considered in the present work is the 
soft repulsive shoulder potential (8). In the works [41. 45| 
it has already been shown that this system demonstrates 
anomalous behavior of diffusion at low temperatures and 
densities corresponding to the region where a competition 
between the characteristic length scales a and ai takes 
place. This competition gives the great complexity of the 
phase diagram of the system [4l|, , so one can expect 
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FIG. 6: Diffusion coefficient of GCM for a set of isotherms, 
(a) - T = 0.04; 0.07 and 0.1; (b) - T = 0.5; 1.0 and 2.0. 



that the thermodynamic quantities, in particular entropy 
which is of the interest of the present study, also have 
a complex behavior in this region of densities. Taking 
into account the complex behavior of both entropy and 
diffusion coefficient it is interesting to check the Rosenfeld 
relation for this system. 

Fig. 11 represents the diffusion coefficient for the soft 
repulsive shoulder system with cti — 1.35 for a set of 
temperatures and densities. One can see that at T = 0.25 
an inflection point in the diffusion coefhcient curve occurs 
which then develops into a loop (T — 0.2). 

Both pair and full excess entropies were considered for 
this potential. Fig. 12 (a) and (b) show the behavior 
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FIG. 7: Radial distribution functions of GCM for a set of 
isotherms at the density p = 0.5 (a) - T = 0.04; 0.07 and 0.1; 
(b) - r = 0.5; 1.0 and 2.0. 



FIG. 8: Excess entropy of GCM for a set of isotherms, (a) 
T = 0.04; 0.07 and 0.1; (b) - T = 0.5; 1.0 and 2.0. 



of the entropies along two isotherms. One can see that 
both at high and low temperature the difference between 
excess entropy and pair contribution to it is rather large. 
This discrepancy is small at low densities, but greatly 
increases at the density about 0.4. Note that this density 
corresponds to a character distance / ~ l/p^^'^ — 1-35, 
that is / ~ CTi . It allows to conclude that the interplay of 
the distances starts at this density and it is this interplay 
which makes the excess entropy and pair excess entropy 
difference to increase rapidly. 

Fig. 13 shows the diffusion coefficient scaling with the 
pair part of the excess entropy. As is seen from the figures 



even at high temperatures the curve is not straight while 
at low temperatures the curve becomes very strange. 
Definitely the exponential relation between the diffusion 
coefficient and pair excess entropy is not valid. 

The diffusion scaling with the full excess entropy is 
shown in the figure 14 (a) and (b). One can see from these 
pictures that the scaling rule works good for the temper- 
atures T = 0.5 and T = 0.4 but aheady for T = 0.35 
the deviation from the linear behavior occurs. This de- 
viation develops more as the temperature decreases. At 
T = 0.25 a self crossing loop occurs. This loop even en- 
larges at lower temperatures. It is worth to note that the 
curve at low temperature T = 0.2 consists of two linear 
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FIG. 9; Excess entropy of GCM for a set of isotherms, (a) 
T = 0.04; 0.07 and 0.1; (b) - T = 0.5; 1.0 and 2.0. 



parts connected by the self crossing loop. It means that 
the Rosenfeld formula is valid in some narrow regions, 
but not in the whole range of densities. 

In order to see the reason for the occurrence of the 
self crossing loop we compare the qualitative behavior 
of the diffusion coefficient and excess entropy. For this 
reasons we measure the densities corresponding to the 
minimum and maximum of the diffusion (/0^j„ and p^ax 
correspondingly) and excess entropy (p^j„ and pfnax) 
the lowest temperature we study T = 0.2. The values we 
obtain are: p£„^ = 0.44, p^^^^ = 0.55, p^„„ = 0.40 and 
Pmax — 0-535. One can see that there is a mismatch in 
the location of the extremal points of these two quantities 



FIG. 10: Reduced diffusion coefficient (formula (1)) of GCM 
for a set of eleven isotherms 



0.20 



0.15- 



0.10 



0.05- 



0.00 



0.05 
0.04 
0.03 
0.02 
< 0.01 




0.3 0.4 0.5 O.E 0.7 0.8 



7=0.50 
7=0.40 
T- 7=0.35 
7=0.30 
•7=0.25 
— 7=0.20 



0.3 0.4 0.5 0.6 0.7 
P 



0.8 



FIG. 11: Diffusion coefficient for soft repulsive shoulder sys- 
tem with (Ti = 1.35 along several isotherms. The inset en- 
larges the isotherms T = 0.5 (squares) and T = 0.25 (circles). 



and therefore there are some regions where the qualita- 
tive behavior of diffusion and excess entropy is opposite. 
Taking into account exponential dependence supposed in 
the Rosenfeld formula one can expect that even small dis- 
crepancies in the qualitative behavior can lead to large 
errors in the scaling relation. 
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FIG. 12: Excess entropy and pair excess entropy for the soft FIG. 13: The diffusion scahng with the pair contribution to 
repulsive shoulder system for (a) T = 0.5 and (b) T = 0.2. the excess entropy at (a) high and (b) low temperatures for 

the soft repulsive shoulder system. 



CONCLUSIONS 

This articles presents a simulation study of the appUca- 
bility of the Rosenfeld entropy scaling to the systems with 
negative curvature and bounded potentials. It was shown 
that the excess entropy scaling can not be applied to such 
systems at low enough temperatures. Interestingly all 
of the systems considered here demonstrate anomalous 
diffusion behavior in some regions of temperatures and 
densities. It makes questionable if the Rosenfeld rela- 
tion is applicable for the systems with diffusion anomaly. 
These results are in contradiction with the results of Refs. 
17, 1^. This contradiction may be attributed to the dif- 



ferences of considered potentials and simulation methods, 
however, this question requires further investigation and 
will be a topic of a subsequent publication. 

One can suppose that the excess entropy scaling is in- 
valid for the systems with negative curvature potentials 
such as repulsive shoulder potential. A possible reason 
for this can be related to the fact that such systems are 
effectively quasibinary j4l|, . As it was mentioned in 
the introduction, the original Rosenfeld idea was based 
on the connection of a liquid under investigation to the 
effective hard spheres liquid. At the same time liquids 
with negative curvature potentials may be approximated 



10 



■1.6n 




■3.2 -2.8 -2.4 -2.0 -1.6 -1.2 



s 

ex 




ui — I — I — I — I — I — I — I — I — I — I — I — |— 
■4.4 -4.0 -3.6 -3.2 -2.8 -2.4 -2.0 



FIG. 14: The diffusion scaling witli the excess entropy at (a) 
high and (b) low temperatures for the soft repulsive shoulder 
system. 



by a mixture of hard spheres of two different sizes. The 
concentration of components of such mixture is pressure 
and temperature dependent. As it was shown in liter- 
ature (see, for example, [!§]) the excess entropy scaling 
holds for binary mixtures too. But in the case of quasi- 
binary mixture since the effective concentration depends 
on the pressure and temperature the behavior becomes 
more complex. This brings to the breakdown of the scal- 
ing rules for this case. 

Obviously, the systems with bounded potentials can 
not be approximated by hard sphere potentials too. It 
seems that this may be the reason of violation of Rosen- 



feld entropy scaling for these systems. 

We thank V. V. Brazhkin and Daan Frenkel for stim- 
ulating discussions. Our special thanks to Prof. Ch. 
Chakravarty who attracted our attention to the prob- 
lems considered here. The work was supported in part 
by the Russian Foundation for Basic Research (Grant No 
08-02-00781). 
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